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ASYMPTOTIC ESTIMATES FOR ROOTS OF THE CUBOID 
CHARACTERISTIC EQUATION IN THE NONLINEAR REGION. 


Ruslan Sharipov 


Abstract. A perfect cuboid is a rectangular parallelepiped. Its edges, its face 
diagonals, and its space diagonal are of integer lengths. None of such cuboids is 
known thus far, though the system of Diophantine equations describing them is easily 
written. The cuboid characteristic equation is a twelfth degree Diophantine equation 
derived from the initial cuboid equations and equivalent to them. In the case of 
the second cuboid conjecture it reduces to a tenth degree equation. This equation 
comprises two parameters. Previously various asymptotics for roots of this equation 
were studied as its parameters tend to infinity either separately or simultaneously 
provided some linear combination of them is preserved finite. In the present paper 
this linear combination is replaced by a certain nonlinear expression. 


1. Introduction. 

Omitting details, which can be found in [1-3], let’s proceed to the tenth degree 
cuboid characteristic equation arising in the case of the second cuboid conjecture: 

Qp,(i) = 0. (1.1) 

The tenth degree polynomial Qpq{t) in (1.1) is given by the explicit formula 

Qpq{t) = + (2 +p2) (3^2 _ 2p2) ^8 (^8 qe+ 

+ - 14p® +P®) (g® - 14p^ g® + dp"' g^+ (1.2) 

+ 10p®g2 +p®)t^ -p®g® (g2 + 2p2) (3^2 - 2 q^) - q^°p^°. 

More details concerning the polynomial (1.2) and its background can be found in 
[4-8] . For the history and various approaches to the problem of perfect cuboids the 
reader is referred to [9-55]. The papers [56-68] constitute a separate stream using 
a symmetry approach to the perfect cuboid equations. This approach is different 
from the approach of the present paper. Therefore we do not consider the results 
of the papers [56-68] below. 

The tenth degree polynomial (1.2) is related to the perfect cuboid problem 
through the following theorem (see Theorem 8.1 in [1] or in [2]). 
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Theorem 1.1. A triple of positive integer numbers p, q, and t satisfying the equa¬ 
tion (1.1) and such that p q are coprime produces a perfect cuboid if and only if 
the following inequalities are fulfilled: 

t>p'^, t>pq, t>q‘^, {p'^+t){pq + t) >2t^. 

The mechanism associating the numbers p, q, t with perfect cuboids was found 
in [4] and [5] . It is described in brief in [1] , [2] , and [3] . We shall not reproduce this 
description in the present paper. 



Fig. 1.1 

Relying on Theorem 1.1, in [1], [2] two forms of asymptotics for the roots of the 
tenth degree polynomial equation (1.1) were studied: 


p = const, 
q = const. 


q -)► +00, 

p —>■ + 00 . 


(1.3) 
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As a result of studying limits (1.3) in paper [2] the following three regions in the 
positive quadrant of the pq - coordinate plane were defined: 

1) linear region given by the linear inequalities 

^<_p, p<59q; (1.4) 

59 

2) nonlinear region given by the nonlinear inequalities 

59q^p, p^9q^; (1.5) 

3) no cuboid region which is the rest of the positive pg- quadrant. 

These regions are schematically shown in Fig. 1.1 above. The linear region (1.4) is 
shown in sky blue. The nonlinear region (1.5) is shown in faded pink. And finally, 
the no cuboid region is shown in gray. 

The linear region was studied in [3] . As a result a narrow strip surrounding the 
bisector line p = q was cut off from the linear region. It is given by the inequalities 

q - —^p, p ^ q + mm(\ (1.6) 

^97 V97 V 74/ ^ ' 

The strip (1.6) was annexed to the no cuboid region. Unfortunately the rest of the 
linear region (1.4) still remains uncertain. Perfect cuboids can potentially be found 
in it, but none of them is actually found. The main goal of the present paper is to 
study the nonlinear region (1.5). 

2. Nonlinear transformation of parameters. 

The upper boundary of the nonlinear region is given by the cubic parabola 
p = 9q^ (see (1.5) and Fig. 1.1 above). For this reason we consider the following 
cubic transformation of the parameters p and q: 

p = Bq^-p, q = q. (2.1) 

Here B is some positive integer number. The transformation (2.1) is invertible: 

p = Bq^-p, q = q. (2.2) 

The transformations (2.1) and the transformation (2.2) map the integer pq-grid 
onto the integer pp-grid and vice versa. 

Note that the curve given by the condition p = const is a cubic parabola. In 
particular, if H = 9 and p = 0, it coincides with the upper boundary of the nonlinear 
region in Fig. I.l. For this reason, instead of (1.3), we set 

p = const, q —>■ +00 (2.3) 

In order to study the limit (2.3) let’s substitute (2.2) into the polynomial (1.2). 
As a result we get another polynomial Qpq{t). This polynomial is given by an 
explicit formula. However, the formula for the polynomial Qpq(t) is rather huge. 
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In the fully expanded form it has 108 summands. For this reason it is placed to the 
ancillary file strategy_formulas_04.txt in a machine-readable form. 

Using the polynomial Qpq{t), one can write an equation similar to (1.1): 

= 0. (2.4) 

It is clear that the equation (2.4) has the same roots as the equation (1.1), though 
they are expressed through different parameters. Like Qpqit), the polynomial 
Qpq{t) in (2.4) is even with respect to t. Along with each root t it has the opposite 
root —t. Therefore we use the condition 

f t > 0 if t is a real root, 

^ (2.5) 

( Im(t) >0 if t is a complex root, 

in order to divide the roots of the equation (2.4) into two groups. We denote 
through ti, ^ 2 , ^ 3 , ^ 4 , ^5 the roots that obey the conditions (2.5). Then te, tr, tg, 
tg, tio are opposite roots of the equation (2.4): 


te — —tl, tr — —t2, ts — — tsi tg — — 14, tiQ — 

3. Parabolic expansions. 

Typically, asymptotic expansions for roots of a polynomial equation look like 
power series (see [69]). By analogy to (2.3) in [1] and according to (2.3), we write 

CX) \ 

l + as g^+cx). (3.1) 

The coefficients Ci in (3.1) should be nonzero: Ci ^ 0. The expansions (3.1) here 
are called parabolic expansions since they occur along cubic parabolas in the plane 
of the original parameters p and q. 

The main tool for studying asymptotic expansions of the form (3.1) for roots of 
a polynomial equation is the Newton polygon of the corresponding polynomial. 

Definition 3.1. For any polynomial of two variables P{t,q) the convex hull of all 
integer nodes (m, r) on the coordinate plane associated with nonzero monomials 
Amr of this polynomial is called the Newton polygon of P{t,q). 

Remark. Note that in our case the polynomial Qpq{t) depends on three variables 
p, q, and t. However, we treat p as a parameter and consider Qpq{t) as a polynomial 
of two variables when applying Definition 3.1 to it. 

The Newton polygon of Qpq{t) is shown in Fig. 3.1 below. Its boundary consists 
of three parts — the upper part, the lower part, and the vertical part. The upper 
part is drawn in green, the lower part is drawn in red. Here are the coefficients of 
monomials associated with the nodes on the upper part of the boundary: 

Ao40 = —.6^°, ^236 = —A432 = — 

(3.2) 

A6 24 = R®, ^8 12 = —2 Aioo = 1- 

Theorem 3.1. The values of exponents ai in the expansion (3.1) for roots of the 


U{p,q) = C,q^ 
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equation (2.4) are determined according to the formula oti = —k, where k stands 
for slopes of segments of the polygonal line being the upper boundary of the Newton 
polygon in Fig. 3.1. 

Theorem 3.1 is a standard fact in Newton polygons application to studying 
bivariate polynomials. Its proof was given in [1] for the sake of reader’s convenience. 



Fig. 3.1 

In our particular case, applying Theorem 3.1, we get 

ai = 2, ai= 4, ai = 6. (3.3) 

The options (3.3) determine the growth rates for the roots of the equation (2.4) as 
q —?► + 00 . They grow as the second, the fourth, and the sixth powers of q. 

The case = 2. This case corresponds to the topmost segment on the upper 
boundary of the Newton polygon in Fig. 3.1. This segment comprises three nodes 
Ao 4 o, ^ 236 j £md ^4 32 - Therefore, substituting the expansion (3.1) with = 2 
into the equation (2.4), we get the following equation for Cp. 


^4 32 + ^2 36 + ^0 40 — 0. 


(3.4) 
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Taking into account (3.2), the equation (3.4) is transformed to 

+ + = 0. (3.5) 

Since i? ^ 0 in (2.1) and (2.2), we can cancel 5^° in (3.5) and write (3.5) as 

Ci^ + eCi^ + l = 0. (3.6) 

The equation (3.6) has four purely complex roots: 

C, = [V2 + l)x, C, = (^-l)i, (3.7) 

a = -(V2 + l)i, C, = -(^2-1)1. (3.8) 


Here i = \/—1. The roots (3.8) are excluded by the condition (2.5). The remain is 
two root (3.7) of multiplicity 1. They yield the following asymptotic expansions: 


ti{p,q) = (^2 + l)ig^ 


ti{P,q) = {V2-l)xq^ 




(3.9) 


The case ai = 4. This case corresponds to the middle segment in the upper 
boundary of the Newton polygon in Fig. 3.1. It comprises two nodes A 432 and 
^6 24 - Therefore, substituting the expansion (3.1) with ai = 4 into the equation 
(2.4), we get the following equation for Ci'. 

^6 24 + ^4 32 = 0. (3.10) 

Applying (3.2), the equation (3.10) is transformed to 

Ci'^ = 0. (3.11) 


Since Ci ^ 0 in (3.1) and since B ^ 0 in (2.1) and (2.2), we can cancel the common 
divisor Ci'^ of two terms in (3.11). As a result (3.11) takes the form 

C^^-B^ = 0. (3.12) 

The equation (3.12) has two real roots, which are simple: 

a =B, Ci = -B. (3.13) 


The second root (3.13) is excluded by the condition (2.5). The remain is one simple 
positive root. It yields the following asymptotic expansion: 

( 00 

l + ^l3isq~" 

S—1 


(3.14) 
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The case Ui = 6. This case corresponds to the lowermost segment in the upper 
boundary of the Newton polygon in Fig. 3.1. It comprises three nodes ^ 524 , 12 , 

and Aiq q- Therefore, substituting the expansion (3.1) with Oi = 6 into the equation 
(2.4), we get the following equation for Ci'. 

AlO 0 + Ag 12 C,® + ^6 24 C,® = 0. (3.15) 

Applying the formulas (3.2), the equation (3.15) is transformed to 

+ = 0. (3.16) 

Since Ci 7 ^ 0 in (3.1), we can cancel the common divisor Ci® of the three terms in 
(3.16). As a result the equation (3.16) takes the form 

Ci'^+B^ =Q. (3.17) 

The equation (3.17) has two real roots of multiplicity 2: 

C^=B^, C,=-B^. (3.18) 


The second root (3.18) is excluded by the condition (2.5). The remain is one simple 
positive root. It yields the following asymptotic expansion: 


U{p,q) 


B^q^(l + J2> 




(3.19) 


The results (3.9), (3.14), (3.19) are summed up in the following theorem. 

Theorem 3.2. For sujficiently large positive values of the parameter q, i. e. for 
q > 9min; the tenth-degree equation (2.4) has five roots of multiplicity one satisfying 
the condition (2.5). Three of them ti, t 2 , and tg are real roots. Their asymptotics 
as q ^ +00 are given by the formulas 


h^B^q^, t2~B'^q^, ts^Bq*. (3.20) 

The rest two roots t^ and tg of the equation (2.4) are complex. Their asymptotics 
as q ^ +00 are given by the formulas 

t 4 ^ {V2+l)iq^, t 5 ^ {V2-l)iq^. (3.21) 


The complex roots (3.21) do not provide perfect cuboids. However, they are 
important for determining the exact number of real roots in asymptotic intervals. 


4. Asymptotic estimates for real roots. 


Acting by analogy to [1], we replace the series expansion of the form (3.1) by 
hnite sum expansions with remainder terms. In the case of the fast growing root 
ti we replace the expansion (3.19) by the following sum: 


ti = B'^ q^+ 2Bq^ 


2Bpf-2q^- 
2pq+ p^ + — - 


20 

52 g2 


+ Ri{p, q,B). 


(4.1) 
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The formula (4.1) is in agreement with (3.19). It means that 


Pii — 0, 

/3i2 = 

2 

5’ 

« _ 2p 

^13 - 

/3i4 

. _ 2p 

/3i6 = 

p2 5 

o 

II 

I> 

/3i8 


<^14 ^ 2 ’ 

20 


in the formula (3.19). Our goal is to derive an estimate of the form 

C 


\R\{PA,B)\ < 


q- 


(4.2) 


for the remainder term in (4.1). In order to get such an estimate we substitute 

t = B^f + 2Bt-2Bpq^-2f-2pq + p^ + ^-^^ + ^. (4.3) 

B B^q^ q-^ 

into the equation (2.4). Immediately after that we perform another substitution 
into the equation obtained by substituting (4.3) into (2.4): 


(4.4) 


z 

Upon two substitutions (4.3) and (4.4) and upon removing denominators the equa¬ 
tion (2.4) is written as a polynomial equation in the new variables c and z: 

l&cB^'^ = SQpB^^ + ip{c,z,p,B). (4.5) 


Here ip{c, z,p, B) is a polynomial of c, z, p, and B with integer coefficients. The 
explicit expression for z,p, B) comprises 1612 monomials. Therefore it is placed 
to the ancillary file strategy_formulas_04.txt in a machine-readable form. 

Let’s recall that H is a positive integer constant. If H ^ 10, then the curve 
p = const goes to infinity outside the nonlinear region (see (2.1), (2.3), and Fig. 1.1). 
Therefore B takes the following finite set of values: 

5 = 1,2, 3,..., 9. (4.6) 


For each B in (4.6) assume that c obeys the condition 

10|p| 


5 = 


< c < 0 if p < 0, 


0 < c < if p > 0. 


(4.7) 


The case p = 0 is exceptional. It should be studied separately. 

Since p is integer and since p 0 in our present case, we have the inequality 


\P\ > 1 - 


(4.8) 


Since we study the asymptotics of the roots ti as q^ -boo, we assume that 

q > 20 ^\. (4.9) 
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Let’s apply (4.9) to (4.4). As a result we derive the inequality \z\ ^ 1/20 
Applying this inequality along with the inequalities (4.7) and (4.8) to the polyno¬ 
mial z,p, B), we derive the following estimate for it: 

|^(c,z,p,S)K72|p|i335. (4.10) 

The inequality (4.10) holds for each value of B in (4.6). 

Now we apply the inequality (4.10) to the equation (4.5). If p < 0, it means 
that the right hand side of the equation (4.5) is a continuous function of c that 
varies from —152 |p| B^^ to —8 |p| B^^ while c runs over the negative interval (4.7). 
As for the left hand side of this equation, it is also a continuous function of c 
that monotonically increases from —160 |p| to 0 while c runs over this interval. 
Therefore the equation (4.5) has at least one root within the negative interval (4.7). 

Similarly, if p > 0, the right hand side of the equation (4.5) varies from 8 |p| B^^ 
to 152 \p \ B^^ while c runs over the positive interval (4.7) and hence the equation 

(4.5) has at least one root within this interval. 

The variable c is related to the original variable t through the formula (4.3). 
Therefore from the above considerations we derive the following inequalities for t: 

B^f + 2Bq'^-2Bpf-2q^-2pq+p^ + ^- 

ID 

+ (4.11) 

-2 5 20 . , . „ 

-2q -2pq+p + — - in the case p < 0, 

B B^ q^ 

B^f + 2Bq'^-2Bpf-2f-2pq+p^ + —- 

B 

20 

-^r^<t<B^f + 2Bq^-2Bpf-2q^- (4.12) 

~2 5 20 lOp . , . „ 

-2pq + P in the case p > 0. 

As a result we have proved the following theorem. 

Theorem 4.1. J/p ^ 0, then for each q ^ 20 -y/j^ and for each value of B in 

(4.6) there is at least one real root of the equation (2.4) satisfying the inequalities 
(4.11) or (4.12) respectively. 

Theorem 4.1 means that we have derived the estimate (4.2) with C = 10 |p|/i3^ 
for the remainder term Ri{p,q,B) in the asymptotic expansion (4.1) for at least 
one root of the equation (2.4). 

The root t 2 in (3.20) is handled in a similar way. The asymptotic expansion 
analogous to the expansion (4.1) for it is written as follows: 

t2 = B^ q^ - 2 B q^ - 2 B pf - 2q^ + 

5 20 

-f 2pq+p'^ ~ -R- u 2 ~2 +^ 2 {p,q,B). 

B B^ 

The formula (4.13) is in agreement with (3.19). It means that 


(4.13) 
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P25 = 


2p 


P26 = 


P 


Pn — 0 




52’ 52 ^ 3 ’ B 

The formula (4.3) in this case is replaced by the following one: 


t = f -2Bq‘^ - -2q^ + 2pq + p^ - -- - 


20 


+ ^. (4.14) 


B B'^q'^ g3 

Upon substituting (4.14) and (4.4) into the equation (2.4) and upon removing 
denominators the equation (2.4) is written as an equation very similar to (4.5): 


16cB 37 = +^(c,z,p,B). 


(4.15) 


Here ip{c, z,p, B) is a polynomial of c, z, p, and B with integer coefficients. The 
explicit expression for ^/’(c, z,p, B) comprises 1612 monomials. Therefore it is placed 
to the ancillary file strategy_formulas_04.txt in a machine-readable form. 

Like in the previous case, for each B in (4.6) assume that c obeys the inequalities 


0 < c < if p < 0, 

±5 ^ 

lOlpI 

- < c < 0 if p > 0. 


(4.16) 


Then assume that q obeys the inequality (4.9). Under these assumptions, taking 
into account (4.8), one can derive an estimate similar to the estimate (4.10): 

|V'(c,z,p,H)K72|p|H35. (4.17) 

Now, writing the inequalities 

B'^ f - 2Bf - 2Bpf - 2q^ + 2pq+p^ “4- 

B 

20 

-^^<t<B^q^-2Bq^-2Bpf-2q^+ ( 4 .I 8 ) 

-2 5 20 lOp . , . ^ 

+ 2pq+p - m the case p < 0 , 

B"^ f - 23^"^ - 2Bpf - 2f + 2pq + p‘^ - 

B 

20 lOp 9 R . O 

-^Bq*-2Bpq^- (4.19) 

2 n 5 20 

-2q +2pq + p - — - ,9 in the case p > 0 , 

B B^ 


and then applying the estimate (4.17) to the equation (4.15) and taking into account 
(4.16), we can easily prove the following theorem. 

Theorem 4.2. //p ^ 0, then for each q ^ 20 and for each value of B in 
(4.6) there is at least one real root of the equation (2.4) satisfying the inequalities 
(4.18) or (4.19) respectively. 
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Theorem 4.2 is similar to Theorem 4.1. It means that we have got the estimate 

l^ 2 (p,g,-B)| <where C = (4.20) 

for the remainder term R 2 {p, q, B) of the asymptotic expansion (4.13). 

The root in (3.20) is somewhat different from ti and ^ 2 - The analog of the 
asymptotic expansions (4.1) and (4.13) for this root is written as 

t^=Bq^-pq+^+R3{p,q,B). (4.21) 

Like in (4.2) and (4.20), in this case we shall derive the inverse cubic estimate 

\R3{P,q,B)\<^ (4.22) 

for the remainder term in (4.21). For this purpose we substitute 

t = Bq^-pq+ — + — 4.23 

B q'^ 

into the equation (2.4). Immediately after that we perform the substitution (4.4) 
into the equation obtained by substituting (4.23) into (2.4). Upon removing de¬ 
nominators, the resulting equation can be written in the following form; 

2cB'^^ =32PB'^^ + f{c,z,p,B). (4.24) 


Here f{c,z,p,B) is a polynomial of c, 2 , p, and B with integer coefficients. The 
explicit expression for f(c,z,p,B) comprises 490 monomials. Therefore it is placed 
to the ancillary file strategy_formulas_04.txt in a machine-readable form. 

For each B in (4.6) assume that c obeys the condition 


32 \P\ 
52 


< c < 0 


0 < c < 


32 Ip I 

52 


if p < 0, 
if p > 0. 


(4.25) 


Again, the case p = 0 is exceptional. It should be studied separately. 
Like in (4.9), assume that q obeys the inequality 


q ^ 7 ^\. (4.26) 

Let’s apply (4.26) to (4.4). As a result we derive the inequality \z\ < l/7|p|“^/^. 
Applying this inequality along with the inequalities (4.25) and (4.8) to the polyno¬ 
mial f{c,z,p,B) in (4.24), we derive the following estimate for it: 

|/(c,z,p,5)K26|p|52i. (4.27) 


The inequality (4.27) holds for each value of B in (4.6). 
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Now we apply the inequality (4.27) to the equation (4.24). If p < 0, it means 
that the right hand side of the equation (4.24) is a continuous function of c that 
varies from —58 \p\ to —6 |p| while c runs over the negative interval (4.25). 
As for the left hand side of this equation, it is also a continuous function of c 
that monotonically increases from —64 |p| B^^ to 0 while c runs over this interval. 
Therefore the equation (4.24) has at least one root within the negative interval 
(4.25). Similarly, if p > 0, the right hand side of the equation (4.24) varies from 
6 IpI B^^ to 58 IpI B^^ while c runs over the positive interval (4.25) and hence the 
equation (4.24) has at least one root within this interval. 

The variable c is related to the original variable t through the formula (4.23). 
Therefore from the above considerations we derive the following inequalities for t: 


Bq"^-pq + ^^ ^ 2^3 <t<Bq ^if p < 0, (4.28) 

Bq^-pq+^ <t < Bq'^ - pq+^ + if p > 0. (4.29) 

B B B‘‘q-^ 

As a result we have proved the following theorem. 

Theorem 4.3. Ifp ^ 0, then for each q ^ 7 \/p| and for each value of B in (4.6) 
there is at least one real root of the equation (2.4) satisfying the inequalities (4.28) 
or (4.29) respectively. 

Theorem 4.3 means that we have derived the estimate (4.22) with C = 32 |p|/i3^ 
for the remainder term R^ifp^q^B) in the asymptotic expansion (4.21) for at least 
one root of the equation (2.4). 


5. Asymptotic estimates for complex roots. 

There are two complex roots t^ and is of the equation (2.4) satisfying the con¬ 
dition (2.5). The leading terms of their asymptotics as 9 —>■ -l-oo are given by the 
formulas (3.21). Specifying the first formula (3.21), we write 

U = {V2 + l)\q^ - +Ri{p,q,B). (5.1) 

Here i = \/—l. Our goal is to derive an estimate of the form 

\R4{p,q,B)\<^ (5.2) 

for the remainder term in (5.1). In order to get such an estimate we substitute 


t = (v^ -I- 1) ig^ 


(10-f7V2)i , ci 


(5.3) 


into the equation (2.4). Immediately after that we perform the substitution (4.4) 
into the equation obtained by substituting (5.3) into (2.4). Upon removing denom¬ 
inators, the resulting equation can be written in the following form: 


30 = -32(10-f 7y2)pH37 + ^(c,z,p,H). 


16ci? 


(5.4) 
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Here r]{c, z,p, B) is a polynomial of c, z, p, and B. The fully expanded expression 
for r](c, z,p, B) comprises 1031 monomials. Therefore it is placed to the ancillary 
file strategy_formulas_04.txt in a machine-readable form. 

Note that the value of the irrational coefficient in (5.4) obeys the inequalities 

636 < 32 (10 + 7V2)^ 636.78 < 640 = 16 • 40. (5.5) 


Therefore for each B in (4.6) assume that c obeys the condition 


0 < c < if P < 0, 

< c < 0 if p > 0. 


80 IpI 


(5.6) 


Like in (4.7), (4.16), and (4.25) above, the case p = 0 is exceptional. It should be 
studied separately. 

Assume that q obeys the following inequality similar to (4.9) and (4.26): 


g ^ 15 ^\. (5.7) 

Let’s apply (5.7) to (4.4). As a result we derive the inequality \z\ ^ 1/15 \p\~^^^. 
Applying this inequality along with the inequalities (5.5) and (4.8) to the polyno¬ 
mial r](c, z,p, B) in (5.4), we derive the following estimate for it: 

|p(c,^,p,H)K512|p|h27. (5.8) 


Let’s compare (5.8) with (5.5) and then apply the inequality (5.8) to the equa¬ 
tion (5.4). If p < 0, it means that the right hand side of the equation (5.4) is a 
continuous function of c that varies from 124 |p| B^'^ to 1152 |p| B^'^ while c runs 
over the positive interval (5.6). As for the left hand side of this equation, it is also a 
continuous function of c that monotonically increases from 0 to 1280 |p| while c 
runs over this interval. Therefore the equation (5.4) has at least one root within the 
positive interval (5.6). Similarly, if p > 0, the right hand side of the equation (5.4) 
varies from —1152 |p| B^^ to —124 |p| B^'^ while c runs over the negative interval 
(5.6) and hence the equation (5.4) has at least one root within this interval. 

The variable c is related to the original variable t through the formula (5.3). 
Therefore from the above considerations we derive the following inequalities for t: 


{V2 + l)q^ - < Imt < (a/2 -|- 1) - 


B^q^ 
10-^772 80p 


(5.9) 


B2q2 B3q5 


in the case p < 0, 


r 10-F 7^2 80p ^ . 

iV2+l)q - <Imt< (V2 + l)g - 

IO + 7V2 


(5.10) 




in the case p > 0. 


The inequalities (5.9) and (5.10) are analogs of the corresponding inequalities in 
the case of the real roots, see (4.11) and (4.12), (4.18) and (4.19), or (4.28) and 
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(4.29). They mean that we have proved the following theorem. 

Theorem 5.1. If p 0, then for each g ^ 15 and for each value of B in 

(4.6) there is at least one purely imaginary root of the equation (2.4) satisfying the 
inequalities (5.9) or (5.10) respectively. 

Theorem 5.1 means that we have derived the estimate (5.2) with C = 80 \p\/B^ 
for the remainder term R^fpjq, B) in the asymptotic expansion (5.1) for at least 
one root of the equation (2.4). 

The complex root t^ is similar to the root t^. For this root we write 

ts = (V2- l)ig^ + +R5{p,q:B). (5.11) 


Here i = V—1- Like in (5.1), our goal is to derive an estimate of the form 

|i?5(p,g,i?)| <-^ (5.12) 

for the remainder term in (5.11). In order to get such an estimate we substitute 


t = {V2-l)\f + 


(10- 7V2)i 

52 g2 


ci 



(5.13) 


into the equation (2.4). Immediately after that we perform the substitution (4.4) 
into the equation obtained by substituting (5.13) into (2.4). Upon removing de¬ 
nominators, the resulting equation can be written in the following form: 


16c 53° = 32(10-7^2)^527 + ((c,z,p, 5). (5.14) 


Here (C(c, z,p, 5) is a polynomial of c, z, p, and 5. The fully expanded expression 
for ({c, z,p, B) comprises 1031 monomials. Therefore it is placed to the ancillary 
file strategy_formulas_04.txt in a machine-readable form. 

Note that the value of the irrational coefficient in (5.14) obeys the inequality 

3 < 32(10-7^2) Ri 3.22 < 4 = 16 • i. (5.15) 

The number in (5.15) is substantially smaller than the number in (5.5). Therefore, 
instead of (5.6) we write the following inequalities: 


IpI 

< c < 0 if P<0, 
0 < c < if p>0. 


(5.16) 


The inequalities (5.16) should be fulfilled for each value of 5 in (4.6). 

Assume that q obeys the following inequality similar to (4.9), (4.26), and (5.7): 


q ^ 3600 


( 5 . 17 ) 






ASYMPTOTIC ESTIMATES FOR ROOTS ... 


15 


The coefficient 3600 in (5.17) is much greater than the coefficient 15 in (5.7). Ap¬ 
plying (5.17) to (4.4), we derive \z\ ^ 1/3600 Applying this inequality along 

with the inequalities (5.15) and (4.8) to the polynomial /'(c, z,p, B) in the equation 
(5.14), we derive the following estimate for it: 

|C(c,z,p,i?)| ^2|p|i?27. (5.18) 

Using the estimate (5.18), we can easily prove that under the assumption (5.17) 
the equation (5.14) has at least one root obeying the corresponding inequalities 
(5.16). The variable c is related to the original variable t through the formula 
(5.13). Therefore from (5.16) we derive the following inequalities for t: 


{V2-l)f + 


10-7^2 


10-7^2 . 


-F 


2B3g5 


< Imt < (\/2 — 1) -F 


(5.19) 




in the case p < 0, 


l)q^ + 2 ^ ^ {V2- l)g^-F 


B^q^ 
10-7^2 p 


(5.20) 




2B3g5 


in the case p > 0. 


As a result we have proved the following theorem. 

Theorem 5.2. //p yF 0, then for each q ^ 3600 y/\^\ and for each value of B in 
(4.6) there is at least one purely imaginary root of the equation (2.4) satisfying the 
inequalities (5.19) or (5.20) respectively. 

Theorem 5.2 means that we have derived the estimate (5.12) with the coefficient 
C = |p|/(2i?3) for the remainder term R^(p,q^B) in the asymptotic expansion 
(5.11) for at least one root of the equation (2.4). 


6. Non-intersection of asymptotic sites. 

In the previous two sections we have used four inequalities (4.9), (4.26), (5.7), 
and (5.17). The inequality (5.17) is the strongest of them. It looks like 

9 ^ 3600 ^. ( 6 . 1 ) 

If the inequality (6.1) is fulfilled, then all of the inequalities (4.9), (4.26), (5.7), 
(5.17) are also fulfilled. Due to these inequalities, applying Theorems 4.1, 4.2, 4.3, 
5.1, and 5.2, we hnd that there are five asymptotic sites, each of which comprises 
at least one root of the equation (2.4). 

The site comprising the root t^ is given by the inequalities (4.28) and (4.29). 
From (6.1) and (4.8) we derive the following inequalities: 


f > 3600^ Ip I 

a3 

fr ^ 3600^ 
\P\ 


q ^ 3600, 

M < ^ 

g3 ^ 36003' 


( 6 . 2 ) 
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From (4.6) we easily derive the following inequalities for B: 
1 ^ B ^ 9, 


1 1 
- ^ ^ 1 . 

9 ^ B ^ 


(6.3) 


Applying (6.2) and (6.3) to the left hand side of (4.28), we derive 


„ _4 __ 16 32p I \ 16 

Bi -p<i + -b + b^s»b, (l-j^)+-g 


32 


^ 3600^ 1 - 


36003 




36003) 


16 32 

^~9~ 36003 


(6.4) 


1.68- 10^'^ > 0. 


36003 

The left hand side of (4.29) is treated similarly: 


i 


> 3600^ 1 - 


3600^ 


1 fi 

— Ri 1.67- lO^"^ > 0. 
9 


(6.5) 


The root t 2 is delimited by the inequalities (4.18) and (4.19). Let’s compare the 
left hand side of (4.18) with the right hand side of (4.28). For their difference, using 
(6.2) and (6.3), we derive the following estimate: 

(^B^q<^-2Bq^-2Bpq^-2f + 2pq+P^-^-^:^)- 


2 p 


2 3 p 

+ 


21 


Pq+ B^f(l ^,2 Bq3 B^q^'B^q^ 

2^\>B^dHi-— _?_?_^ 

V ^ ^ V 36002 30003 30004 30095 

) ft! 0.99 - B ^ f ^ 0.99 • 3600® ft 2.18 • 10^^ > 0. 


( 6 . 6 ) 


B^q& B'^q^, 
21 20 


36006 30008 

The difference of the left hand side of (4.19) and the right hand side of (4.29) is 
treated similarly. For this difference we have 


{B'^q^ -2Bq^ -2Bpf - 2 f + 2pq + p^ - ^ - 

21 


20 


52^2 


lOp 
B'^ q^ 
2 


3600^ 


20 


42 


(6.7) 


36003 30004 30005 30006 30008 30000 / 

ft 0.99 • g® > 0.99 • 3600® ft 2.17 • 10^^ > 0. 

The root ti is delimited by the inequalities (4.11) and (4.12). Let’s consider the 
difference of the left hand side of (4.11) and the right hand side of (4.18): 


(B^q^ + 2Bq^-2Bpf-2e-2pq+p^ + ^-^^ + 

+ - (i?"9®-2Sg"-2i?pg®-2g2 + 2pg + p2_ 
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20 


10 p 


B B‘^q^ B^q^ 


) =4Bq^-4pq+-^ + 


10 20 p 


B'^q^' 

Applying (6.2) and (6.3) to the above expression, we get the following estimate: 

p 5 


10 20p 

ABq —4:pq + — + 


B ' B^q^ 

B^q^J ^ ^ V 360 


Bq^ 2B2g4 

5 5 \ 


36003 2-36004 3600^ 


( 6 . 8 ) 


4.00 -Bq^^^ 4.00 • 3600"^ « 6.71 • 10^^ > 0. 


The difference of the left hand side of (4.12) and the right hand side of (4.19) is 
treated similarly. For this difference we have 

(i ?2 56 + 2594 -2i?p93_2g2_2pg+p2 + |_^^)_ 


-(52 56_2B54_25p93_292 + 2p9 + p2_|_^^) 


= 4:Bq^ - ipq+ — ^ ABq'^fl - 

y py-r ^ 


36003 2•36004 

4.00 • 5^4 ^ 4.00 • 36004 _ g . jq14 ^ Q 


) 


(6.9) 


The inequalities (6.4), (6.5), (6.6), (6.7), (6.8), and (6.9) mean that if the inequality 
(6.1) is fulfilled, then the asymptotic sites for the roots ti, t 2 , ts do not intersect 
with each other and are located within the positive half-line of the real axis. 

Now let’s proceed to the purely imaginary roots and The root is delimited 
by the inequalities (5.19) and (5.20). Applying (6.2) and (6.3) to the left hand side 
of (5.19), we derive the following estimate for it: 


(72-1)9^-^ 


10-7^ 


X 1-F 


5292 25395 

3^2-4 , piy/2 + 1) 

B'^q'^ 

p(y2 + l) 


25397 




= [V2-l)fx 

3^2-4 


(72 - 1) f (1 - 


36004 


( 6 . 10 ) 


) ft! 0.41 • f ^ 0.41 • 3600^ ft 5.31 • 10® > 0. 


2-36007 

The left hand side of (5.20) is treated similarly: 

10-772 


(72-1)9^ + 
s(72-1)9"( 


5292 

372-4 

' 36004 


(72-1)9'( 


1-F 


372-4 


) 




5294 

) > 0.41 -3600^ ft 5.31 - 10® > 0. 


( 6 . 11 ) 


The root ti is delimited by the inequalities (5.9) and (5.10). Let’s compare the 
left hand side of (5.9) with the right hand side of (5.19). Here is their difference: 


2 +1) r - 


2 10-F772\ 


5292 


2 - 1 ) 9 " + 


2 , 10-772 \ 20 


52 92 


-) =29'"- 


5292 - 
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Applying (6.2) and (6.3) to the above expression, we get the following estimate: 


2^- 


20 


^ \ B^q^J ^ \ 3600 ^^ 


10 




( 6 . 12 ) 


1.99 • > 1.99 • 3600^ Ri 2.59 • 10^ > 0. 


The difference of the left hand side of (5.10) and the right hand side of (5.20) is 
treated similarly. For this difference we have 


2 10 + 7^2 


2 + 1)9 Q2q2 

10-7^2 p 


80 p 


B^q^ 

2 B^q^) = 29 "- 

10 81p 


20 


2-l)g" + 

81 p 


B^q* 2B^q 


1 - 


B^q^ 

= 29 ^( 1 - 

« 1.99 • f ^ 1.99 • 3600^ 2.59 • 10^ > 0. 


B"^ ip B^ p 

) \ „ -2 / 10 


3600^ 


(6.13) 


The inequalities (6.10), (6.11), (6.12), and (6.13) mean that if the inequality (6.1) is 
fulfilled, then the asymptotic sites for the roots Ia and do not intersect with each 
other and are located within the positive half-line of the imaginary axis. Summa¬ 
rizing this result with the above result for the real roots ti, t 2 , P we can formulate 
the following two theorems. 

Theorem 6.1. // p < 0 and q ^ 3600 y/p|, then the equation (2.4) has five simple 
roots obeying the condition (2.5). Three of them are real and positive. These three 
positive real roots are located within three disjoint asymptotic sites given by the 
inequalities (4.11), (4.18), and (4.28) respectively. 

Theorem 6.2. If p > 0 and q ^ 3600 \/p|, then the equation (2.4) has five simple 
roots obeying the condition (2.5). Three of them are real and positive. These three 
positive real roots are located within three disjoint asymptotic sites given by the 
inequalities (4.12), (4.19), and (4.29) respectively. 


7. Integer points of asymptotic sites. 

Let’s consider the inequalities (4.11) and (4.12) delimiting the root ti of the 
equation (2.4). Most of the terms in them are integer, except for following ones: 


5 20 


lOp 
B'^ p 


(7.1) 


The first term (7.1) is optionally non-integer. If the inequalities (6.2) and (6.3) are 
fulfilled, then the other two terms (7.1) are certainly non-integer. 

Let’s begin with the case where 5/i? is not integer. In this case due to (6.3) it 
is separated from the nearest integer number by a distance not less than 1/9: 




(7.2) 
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The other two terms in (7.1) are substantially smaller. From (6.2) and (6.3) we get 


20 




< 


20 

36002’ 


< 


10 


lOp 

B2g3 - 30003- 


Combining (7.2) and (7.3), we obtain the following inequality: 


4± 


20 


IQp 

± .o — n 


B B'^q^ 


1 

^ 9 


20 

36002 


10 


1 

> TTT- 


36003 10 


(7.3) 


(7.4) 


In the case where 5/i? is integer, i. e. where 5 = 5, we treat the inequalities 
(4.11) and (4.12) separately. If p < 0, from (6.2) and (6.3) we derive 


-1 < -- 


20 


10 


^ -- 


20 


lOp 

36002 30003 ^ ^2^2' 52^3 

If p > 0 we need an additional condition for p and g: 


< - 


20 


52 g 2 


< 0 . 


(7.5) 


2g > IpI 


(7.6) 


Provided the condition (7.6) is fulfilled, from (6.2) and (6.3) we derive 
20 ^ 20 ^ 20 IQp 

^ “ 3 ^ ^ ^ ^ ^ 

If the condition (7.6) is not fulfilled, then from (6.2) and (6.3) we derive 
20 ^ 20 ^ 20 lOp ^20 ^ 

^ ~36002 ^ 2 q 2 < “52^2 + 52^3 < 30002 + 30003 < 


(7.7) 


(7.8) 


Theorem 7.1. If p ^ 0, then for each q ^ 3600 y/pj and for each B b in (4.6) 
the asymptotic site given by the inequalities (4.11) and (4.12) has no integer points. 

Theorem 7.1 is proved by applying (7.4) to (4.11) and (4.12). 

Theorem 7.2. If p < 0 and 5 = 5, then for each q ^ 3600 y/p| the asymptotic 
site given by the inequalities (4.11) has no integer points. 

Theorem 7.2 is proved by applying (7.5) to (4.11). 

Theorem 7.3. If p > 0 and 5 = 5, then for each q such that q ^ 3600 \/p| and 
2q> p the asymptotic site given by the inequalities (4.12) has no integer points. 

Theorem 7.3 is proved by applying (7.7) to (4.12). 

Theorem 7.4. If p ^ 0 and 5 = 5, then for q ^ 3600 \/pj the asymptotic site 
given by the inequalities (4.11) and (4.12) has at most one integer point given by 
the formula f = 25g® + lOg^ — lOpg^ — 2q^ — 2pg + p2 + 1, where p > 0. 

Theorem 7.4 is proved by applying (7.8) to (4.12). 

Note that the inequalities (4.18) and (4.19) are quite similar to (4.11) and (4.12). 
Applying the estimates (7.4), (7.5), (7.7), and (7.8) to them we can prove the 
following four theorems for the corresponding asymptotic site. 
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Theorem 7.5. If p 7 ^ 0, then for each q ^ 3600 and for each B ^ 5 in (4.6) 
the asymptotic site given by the inequalities (4.18) and (4.19) has no integer points. 

Theorem 7.6. If p > 0 and B = 5, then for each q ^ 3600 \/p| the asymptotic 
site given by the inequalities (4.19) has no integer points. 

Theorem 7.7. If p < 0 and B = 5, then for each q such that q ^ 3600 \/p| and 
2q> IpI the asymptotic site given by the inequalities (4.18) has no integer points. 

Theorem 7.8. // p 7 I 0 and B = 5, then for q ^ 3600 \/[^ the asymptotic site 
given by the inequalities (4.18) and (4.19) has at most one integer point given by 
the formula t = 25 q^ — 10 q'^ — lOpq^ — 2 q^ + 2pq + p^ — 1, where p < 0. 

Let’s proceed to the inequalities (4.28) and (4.29) and let’s begin with the case 
where 16/5 is not integer. In this case due to (6.3) the number 16/5 is separated 
from the nearest integer number by a distance not less than 1/9: 


16 

'B 


— n 




(7.9) 


From (6.2) and (6.3) for the 


other fractional term in (4.28) and (4.29) we derive 


32p ^ 32 

52^3 -- 36003 


Combining (7.9) and (7.10), we obtain the following inequality: 


(7.10) 


16 , 32p 1 32 1 

5 52 g 3 ^ 9 36003 ^10 


(7.11) 


The case where 16/5 is integer is more simple. In this case we write the inequality 
(7.10) in the following slightly modified form: 


0 < 


32 p 

52 g 3 




32 


3600' 


< 1 . 


(7.12) 


In both cases, applying either (7.11) or (7.12) to the inequalities (4.28) and (4.29), 
we can prove the following theorem. 

Theorem 7.9. If p ^ 0, then for each q ^ 3600 \/Pl and for each 5 in (4.6) the 
asymptotic site given by the inequalities (4.28) and (4.29) has no integer points. 


The exceptional case. 

As we noted above the case p = 0 is exceptional, see Theorems 7.1 through 7.9. 
This case should be studied separately. Let’s substitute p = 0 into the equation 
(2.4). The resulting equation can be written explicitly: 

t^° + (6g^- 2 5^g^^ - f B^)f + {B^f^ + lOq^^B^+Aq^^B^- 

- Uq ^° B ^ + f ) + (14 ^4 _ 4 -24 ^6 _ ~16 ^2 _ -32 ^ 10 _ 

- 10^8) ^4 (2 5® - 6g®® B^° + g®^ B®) f - q'^° B^° = 0. 


( 8 . 1 ) 
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Theorem 8.1. For each B in (4.6) the polynomial in the left hand side of the 
equation (8.1) is irreducible in the ring Z[<:]. 

Theorem 8.1 is proved by means of direct computations. It means that for each 
B in (4.6) the equation (8.1) has no integer roots. 


9. Application to the cuboid problem. 


Theorems 7.1 through 7.9 are based on the inequality (6.1), where p ^ 0- The¬ 
orem 8.1 means that we can omit the condition p 0. Transforming (6.1) back to 
the initial variables p and q with the use of (2.1), we get the inequalities 


Bq^ - 


3600^ 


ifpi^Bq^ + 


3600- 


(9.1) 


From Theorems 6.2, 7.1, 7.5, and 7.9 we derive the following result. 

Theorem 9.1. For each B ^ 5 in (4.6) if the inequalities (9.1) are fulfilled, then 
the tenth degree cuboid characteristic equation (1.1) produces no perfect cuboids. 

The case B = 5 is special. In this case we have the additional condition (7.6). 
Upon transforming (7.6) back to the initial variables p and q it looks like 


Bq^-2q<p<Bq^+ 2 q. 


(9.2) 


From Theorems 6.2, 7.2, 7.3, 7.6, 7.7, and 7.9 we derive the following result. 

Theorem 9.2. For B = 5 if the inequalities (9.1) and (9.2) are fulfilled, then the 
tenth degree cuboid characteristic equation (1.1) produces no perfect cuboids. 

The inequalities (9.2) do not follow from (9.1). They become very restrictive for 
large g as g —>■ -|-oo. Theorems 7.4 and 7.8 can be applied in order to remove this 
restriction. However, they do not change the state of affairs in general. Therefore 
this step is left for one of the future papers. 


10. Conclusions. 

The main result of this paper is presented by Theorems 9.1 and 9.2. Theorems 9.1 
and 9.2 shrink the nonlinear region on the pg-coordinate plane by cutting off nine 
narrow strips surrounding nine cubic parabolas 

p = B q^, where H = 1, 2, ... , 9. 

These strips are outlined by the inequalities (9.1) and (9.2). They are annexed to 
the no cuboid region (see Fig. 1.1) thus reducing the area where perfect cuboids 
are still potentially possible^. 
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